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Abstract 

We present the renormalization scheme used in and the characteristic features 
of GRACE/SUSY-loop, the package of the program for the automatic calcula- 
tion of the MSSM processes including one-loop order corrections. The two- 
body and three-body decay widths of charginos in one-loop order evaluated 
by GRACE/SUSY-loop are shown. 

1 Introduction 

The concept of supersymmetry (SUSY) [H [2] is considered the most promis- 
ing extension of the Standard Model (SM) [3] of particle physics. Among 
the super symmetric theories, the minimal supersymmetric standard model 
(MSSM) is the most elaborated and well studied framework of SUSY. The 
existence of many new supersymmetric particles in any SUSY model makes 
it very complicated and tedious to compute even a simple two-body decay 
width exactly. To overcome this problem, the Minami-tateya group of KEK 
has constructed a computational system GRACE/SUSY [H [5] which automati- 
cally creates, for a given process, all the Feynman diagrams and compute the 
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Feynman amplitudes and subsequently the cross section or the decay width 
itself at tree level in the MSSM. 

With the increase of more precise experimental data, in particular, in 
future colliders, it becomes apparent that the inclusion of at least one-loop 
radiative corrections is necessary for calculations of cross sections of SUSY 
processes. For this purpose, the SPA project has established the conven- 
tion of the SUSY parameters and loop calculations |6]. We have, therefore, 
extended GRACE/SUSY in order to incorporate radiative corrections in one- 
loop order in the system. The new system, called GRACE/SUSY-loop, is con- 
structed based on the same philosophy as for GRACE-loop [7j developed by 
the Minami-tateya group, which is the automatic computation system for 
the SM processes including one-loop corrections. Some of the results com- 
puted by this system have already been presented in several workshops and 
symposiums [5| [TU]. 

In this paper we present the result of the thorough investigation of the 
two-body and three-body decays of charginos in one-loop order in the MSSM 
computed by the automatic computing system GRACE/SUSY-loop. We also 
show the cross section of the chargino pair production in e^e~ annihilation 
and the subsequent decays by combining the production cross section and 
the decay rates. 

The paper is organized as follows. In section two, the renormaliza- 
tion scheme used in GRACE/SUSY-loop is explained. The features of the 
GRACE/SUSY-loop system are briefly given in section three. The numerical 
results are presented in section four, and several comments on the numerical 
results are given in section five. 

2 Renormalization scheme 

In this section we explain briefly the renormalization scheme adopted in 
GRACE/SUSY-loop. Our approach is a straightforward extension of the on- 
shell renormalization in the SM used in GRACE |7j. 

The Lagrangian of the MSSM has been given in pTl [T2| [T3] . In terms of 
superfields, it is given as ( see [TT] for detail) 
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+ J dWe^iexpig'Y^V ~ g,\"*Vn^u 

+ J d^ed^e^leMg'YdV - 9s>^'"*K")^d 

Vi J 
V2 J 

+^ fd'e^n,^^^d + h.c. 

Vi J 

-fjL J d^6l$Hl*H2 + h.C. 

+^soft 

+^^9/ + ^ghost, (2.1) 

where W, W, Wg are superfield strengths corresponding to the SU(2)x,, U(l) 
and SU(3)c gauge-superfields, V, V, Vg, respectively. The component fields 
belonging to each superfield are given as follows 

V = (Wf^X), 5?7(2)l gauge boson 

V = {B^,\), f/(l) gauge boson 
Vs = {g",g'^), 5[/(3)c gauge boson 



Ha 
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= {dR,d*R), (2.2) 

The SU(2) L doublet gauge-bosons and gauginos are denoted by W^, A", while 
the singlet gauge-boson and gaugino arc denoted by 5^ and A. The Higgs 
and higgsino doublets are denoted by Hj and Hj with i — 1,2, respectively. 
The soft SUSY breaking terms are expressed as 

Csoft = -ImiAA - iMsA^A'^ - ^Ms^r + ^-c. 

-m^HtHi - m^H^Hs - (m^^HiH^ + h.c.) - ^ rh}JUm 

fm 

^'"M„H2A(g^)^(K«) + ^^A,ll^A{qL)A{dR) + h.c. 



V2 Vi 

+::/^A,H^A{h)A{eR) + h.c. (2.3) 

where the sum in the sfermion mass terms runs for = and fji. The 
sign convention of our ^/ is opposite to the convention used by others. 

Sfermion mass eigenstates are denoted by /j with i — 1,2 which are the 
mixture of the left-handed (/l) and the right-handed sfermions (/r). The 
sfermion mass matrix is diagonalized as, 

r cos% sin^;^ (ml mJ^J r cos Of -sinOf^ 
[—sinOf cos6f ) "^/*^ "^/^ J Isin^/ cos6f J 



where 



mjr 

ji 

m| 

72 



(2.4) 



^ +m5 + M| cos 2/3(^3/ -Q/s^), 

= rh'j^ + m'j + Ml cos 2 PQfs'^, (2.5) 

^2 ^ ( -niudi cot P + A^), f^u 
fi-R )—mf{iJ,tanf3 + Af), f — d,e 
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Their masses and the mixing angles satisfy the relations originating from 
the S\J{2)l conditions on their left-handed soft SUSY-breaking mass terms, 
m? = rrij , m? = m? etc. For the third generation of sfermions, for 
example, 

cos^ Otml^ + sin^ 6'tm| — ml = cos^ ^bmf^ + sin^ ^fe^^fj ~ + ^'^^ 2/3, 

m?^ = cos^ Orml^ + sin^ Orml^ — ml + cos 2/5. 

(2.6) 

Regarding now the Lagrangian fl2.1l) and fl2.3p as a bare Lagrangian and 
renormailzing all the quantities, we separate the renormalized Lagrangian 
and its counterterms. The renormalization constants are introduced as fol- 
lows. 

[Standard Model sector] 



gauge bosons Wf,o = Z\l^W^, 

r> '7I/2 75 

_ 7I/2 

gauge couphngs = ZgZ^^'^g, 

I ry ry — 'i/2 

9o — ^a'^B 9 
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fermions ^! jl^ = ^''^'^ jl, / = m, (i, ■ ■ ■ , i/e, e. 



^/i?o = Zf^^'^^ifR, / = M,ci,- ■■,e,- ■■ , 
m/o = mf + 6mf, f = u,d, ■ ■ ■ ,e, ■ ■ ■ . (2.7) 



[SUSY sector] 



Higgs bosons 



sfermions 



H 



iO 



^2)0 



71/2TT 



Vi - 6vi 



Z-^^\ml^bm]) 



i = l,2 

z = l,2^ 
i = 1,2 



mi 
1/2 



/i/i /1/2 
^1/2 ^1/2 

/2/1 /2/2 



(m|)o = 



^v, z/j, i = e,fi,T , 

m| + 6m\ , f = u,d, ■ ■ ■ ,e, ■ ■ ■ , 

Ji Ji 

z = l,2, 



Inos 




Of + 89 f, f = u,d, - ■ ■ ,e, - ■ ■ , 




H 



Mio 
M30 



Ml + 5Mi, 
M2 + 8M2, 

M3 + 5M3, (mp = mg + 8mg) , (2.8) 



where in the Higgs boson part, 



ml = ml + l/i 



2 



« = 



1,2 . 



(2.9) 



Summing up, we have (3 + TNg) wavefunction renormahzation constants, 
{?> + ?>Ng) mass counterterms in the non-SUSY sector, and (7+ 13A^g) wave- 
function renormahzation constants and (9+10A'"g') mass, vacuum-expectation- 
value and mixing-angle counterterms in the SUSY sector. 

Several comments are in order. The relation (12. 4p ~ (12. 6p are originally 
satisfied by bare quantities, but they are also valid among the renormalized 
quantities, although the renormalized masses are not necessarily equal to the 
pole masses. In a similar way, (12.91) can be understood as the relation among 
the renormalized quantities. Note that at this state all the renormalized 
parameters in the Lagrangian are simply parameters of the model, and only 
upon imposing the renormahzation conditions to be specified in this section, 
they are related to the physical quantities. 

In the Higgs and higgsino sectors, the wavefunction renormahzation con- 
stants are introduced to each unmixed bare doublet state. The mixing angles 
in the Higgs, chargino and neutralino sectors are defined as the angles which 
diagonalize the renormalized mass matrices. Therefore, there appear no bare 
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mixing angles or counterterms for the mixing angles of charginos and neu- 
tralinos in our scheme. See |14j . 

In place of 6m1, 6m\ and 5m\2i we use the mass counterterm of the CP 
odd Higgs particle, 5M\, and two counterterms of the tadpole interactions, 
5Ti and 572, which are given by the linear combination of 5m\, 5m\ and 

'■121 



(5m? 



5M\ = sl6ml + cl6ml-2cpsp6ml2 (2.10) 



- sp) i^Zz + 6Zh, + 6Zh, - 2 2 — + 2 2 — 
[5m^ — mlSvi + V26m\2 — m] 



+§(5^?' + 5g"){vl - vl)v, + + g"){{vl - 3vl)6v, + 2v,V26v2} 
+W + 9"){i^v! - v,vI)5Zh, - v,vI5Zh,}, (2.11) 

o 

ST2 = (1^2) in 6T1. (2.12) 

We introduce the gauge-fixing Lagrangian in terms of the renormalized 
fields as we have done in GRACE- loop for SM [7j. No renormalization con- 
stants are introduced for gauge fixing constants and no counterterm La- 
grangian corresponding to the gauge-fixing Lagrangian appears. The renor- 
malization of the ghost fields is not necessary in one-loop order. 



We use the on-shell renormalization scheme. The renormalization con- 
ditions employed in GRACE/SUSY-loop are the following set of conditions. 
Using these conditions, we can express all the renormalization constants and 
counterterms in terms of the linear combination of the two-point functions 
evaluated at some specific renormalization points, 
gauge sector 

We use the on-mass-shell condition for W, Z and photon. In addition, we 
require that the residue of the photon propagator at the pole position is one, 
and that the A^-Z^j, transition vanishes for the on-shell photon. 
SM fermions 

We use the same renormalization conditions adopted in GRACE-loop. Namely, 
we require the on-mass-shell condition and the residue condition that the 
residue of the fermion propagator at the pole is one. 
Higgs sector 

We impose the on-mass-shell and the residue condition for CP odd Higgs, 
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A^, and the decoupling of and A'^ on the mass-shell of A^. For the CP 
even Higgs, we impose the on-mass-shell condition for the heavier Higgs, H^. 
The above three conditions together with one of the conditions imposed in 
the gauge sector determine four renormalization constants 5 Hi, 5H2, ^ and 

Note that we have not adopted the often erroneously used renormalization 
condition, ^ = since this condition violates the gauge invariance |15j . 
tadpole terms 

Identical to GRACE-loop, we require that the tadpole terms in the renor- 
malized Lagrangian vanish by itself and the tadpole counterterms cancel the 
one-loop tadpole contributions, 
chargino sector 

We impose the on-mass-shell condition on both xf and xt- addition, we 
impose the residue condition on x^^. 
neutralino sector 

We impose the on-mass-shell and the residue condition on the lightest neu- 
tralino, Xi- 
sfermion sector 

We impose the on-mass-shell condition and the residue condition on all the 
seven sfermions in each generation. In addition, we impose that there is no 
induced mixing between physical /i and f2- The counterterm for the slepton 
mixing angle is determined by the SU{2)l relation (12. 6p upon introducing 
the renormalization constants (12. 7p and (12.81) . In the squark sector, there are 
two counterterms 66^ and 66d for each generation. We fix 66u [IS] by 

2 - mi^ 

while 66d is fixed by the SI]{2)l relation (12.61) upon introducing the renor- 
malization constants (12. 7p and (12.80 . 
charge renormalization 

The charge (electromagnetic coupling constant) is defined as in the standard 
model [l7j. 
QCD sector 

We impose the on-mass-shell condition and the residue condition for gluons 
and gluinos. The counterterm SZg^ is determined by the minimal subtraction 
with dimensional reduction (DR). 
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The explicit expression of the renormahzation constants in terms of two- 
point functions is given in Appendix A. 

A couple of comments are worthwhile at this stage. In our renormal- 
ization scheme, the pole mass of h^, H^, X3 ^ind X4 is different from 
its Born value. Therefore, the mixing angles of the Higgs sector and the 
neutralino sector are not directly related to the pole masses. They are con- 
sidered as effective parameters. This means, in particular, that a quan- 
tity "tan/?" which includes one-loop corrections is not defined and not used 
in GRACE/SUSY-loop. The relation between the parameter tan/? used in 
GRACE/SUSY-loop and the experimentally observed "tan/?" which includes 
higher order corrections depends on how "tan/?" is actually defined. 

In addition to the wavefunction renormahzation constants introduced in 
fl2.7l) and fl2.8p . we need to introduce the ultraviolet finite external wavefunc- 
tion renormahzation constant for each particle for which the residue 
condition is not imposed on its propagator, namely for W^, H^, h^, H^, 
X21 X25 Xs! X% The expression of 5Z^^* in terms of two-point functions is 
given in Appendix B. 

Since the coupling constant of the soft SUSY-breaking Yukawa interac- 
tion among Higgs and sfermions, Af, appears always in combination with 
m/, we use in GRACE/SUSY-loop mfAf and its counterterm 6{mfAf) as in- 
dependent variables, which helps avoid partly the numerical instability of the 
counterterm at large tan /?. Explicitly 

S{mfAf) = - (5 mj^ — 5 J sm29 f + 59 f{m'^j:^ — mj J cos 29 f 

_ f (5(mt/i cot/?) f = u,c,t 

1 5(mj/i tan/?) t = d, s,b,e, ^ 

since according to (12.41) Af is related to the sfermion masses and m^^ as 
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m,A, = cos0,sin^,(m|-m5-)-{;;;^J^^'°|;^ / = 6, e, r " ^^-IS) 

The system can easily accommodate different renormahzation schemes 
by re-expressing the renormahzation constants and the mass counterterms in 
terms of different linear combinations of the two-point functions. 

A brief comparison of our scheme with other earlier studies is worthwhile 
to clarify the difference and also the possible scheme dependence in different 
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schemes. There are plenty of papers on the renormalization schemes of the 
MSSM, and it is beyond our scope to compare all of them. The earlier works 
on the MSSM renormalization (see for example, [TSl [TU ITH]) are naturally 
concerned with radiative corrections in the Higgs sector. Later, the study 
of the renormalization of the other sector, gaugino and higgsino sector and 
sfermion sector followed ( see for example, [20l|2ll[22l|23]). In[18l|19l[Tll 
[22| [21] the on-shell renormalization is used, while the so-called DR is used in 
[21j . Since the number of the particles exceeds the number of the parameters 
appearing in MSSM, one cannot impose the on-shell conditions on all the 
particles. Therefore, there can be many variants even among the on-shell 
scheme, depending on which particles are put on- mass-shell. 

The renormalization scheme we adopt in GRACE/SUSY-loop is close to the 
scheme given by [H], but there are a couple of differences. In [H], different 
from our scheme, the gauge-fixing Lagrangian is introduced in terms of bare 
fields. Therefore, in their scheme, they need extra renormalization conditions 
to fix the renormalization constants for gauge-parameters. Another difference 
of our scheme from the others lies in the renormalization condition imposed 
on 5vi. The condition used in [IS], [HI [22] 5vi/vi = 6V2/V2 is not preferable on 
the ground of potentially violating the Ward identity. We use the on-mass- 
shell condition for in place of the condition on the vacuum expectation 
value, which leads to our expression of 6 tan f3 which is also different from 
others. 



3 Features of GRACE/SUSY-loop 

We present in this section several features of the current version of the system 
GRACE/SUSY-loop. Some of them have been given in 

In order to check and detect possible errors in the system we have used the 
technique of the non-linear gauge (NLG), by adding the SUSY interactions 
to the gauge fixing functions of the SM [T7j. Explicitly, we used the following 
gauge fixing function^ 

= -^\FwA'-7^iFzr-^{F,r, (3.1) 
c,w ^Kz ^K-t 



^In principle, we can also add non-linear sfermion interactions. We have not attempted 
this extension, since due to the mixing the resultant Lagrangian becomes too lengthy and 
cumbersome. 
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Fw± = id^±ieaA^±tgcwf3Zf,)W^^ 

±i^w^{v + 6HH'^ + Shh° ±iliG^)G^, (3.2) 

Fz = d^Z^ + ^ZYi^ + ~eHH' + ihh')G', (3.3) 

= d,A>', (3.4) 

where 

V = ^vl + vl . (3.5) 
The gauge fixing functions contain now seven free parameters, 

^H, Sh, K, in, h- (3.6) 

Since physical results are independent on the NLG parameters, they must 
vanish in the sum of the Feynman amplitudes for physical processes. The 
test using the NLG parameters provides us with a more powerful tool for the 
check of the system than the test using the linear gauge parameters, because 
each NLG parameter is concerned with many kinds of amplitudes which are 
not in the same gauge- independent sub-set in the linear gauge [TO] . 

In the actual computation of decay widths and production cross sections 
in GRACE/SUSY-loop, we use the 't Hooft- Feynman gauge with = 1. For 
the consistency check of the computation, we use the vanishing of the ul- 
traviolet divergences and the infrared singularities in the sum of the loop 
and the soft photon/gluon radiation diagrams, as well as the stability in the 
sum of the soft photon/gluon radiation and the hard photon/gluon radiation 
diagrams against the change of the photon/gluon energy cut-off. 

The input parameters of GRACE/SUSY-loop are 

e, gs, Mw, Mz, M^o, tan/3, /x. Mi, M2, M3, 

0u,dd,0e, ■ ■ ■ , (3.7) 

Using (12.61) . we fix the remaining two masses of the sfermions in each gen- 
eration. Note that the width of the unstable particles is neglected in the 
computation of the amplitudes. 

The coupling constants ^/ of the soft SUSY-breaking Yukawa interac- 
tion among Higgs and sfermions are not our independent input parame- 
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tan/3 




Ml 


M2 




Mao 


10.00 


399.31 


100.12 


197.52 


610 


424.9 

















cos 9u 


cos 6d 


cos 9e 


545.67 


563.44 


545.50 


569.03 


125.50 


190.14 


172.70 


2.4x10-^ 


0.011 


1.1x10-* 
















cos 9c 


cos 6 s 


COS 9fj_ 


545.66 


563.45 


545.52 


568.97 


125.43 


190.16 


172.69 


0.063 


0.018 


0.023 
















cos 9t 


cos 9b 


cos 9r 


368.53 


583.79 


450.12 


544.38 


107.71 


195.08 


170.63 


0.722 


0.967 


0.314 



Table 1: The value of the MSSM input parameters for set (A) 



ters, since they are expressed in terms of sfermion masses and the mix- 
ing angles as fl2.15p . We don't use the GUT relations Mi = | tan^ 6'vi^M2, 

Ms = |t sin^ 6wM2 at our energy scale. 

We have chosen the input values of the parameters in such a way that 
they reproduce the the SPAl pole mass values [6] as close as possible. Since 
the input values of the SPAl parameters |6] are defined by the DR scheme, 
it is not possible to reproduce exactly the same values as those proposed by 
SPAl. Note that the lightest neutralino (x?) is the lightest SUSY-particle 
(LSP) both in SPAl and in our parameter choice. 

In the numerical calculation we adopt two numerical sets (A) and (B). 
which are given in Table 1 and 2, respectively. For the fermion masses of the 
SM and gauge boson masses, we used the following values (in unit of GeV) 
both for sets (A) and (B), 



Miy = 80.35, Mz = 91.1876, 
0.51099906 X 10"^ 



m 



105.658389 x 10 



1.5, 



rrie = 
rriu -- 

nid = C)».U X lu rris = 'dZV x lu 
while for the strong coupling constant, we used 

a, = 0.12 . 



-3 



1.7771, 



58.0 X 10-^ 
58.0 X 10~^ 



= 0, 
rric 



(3.8) 



mt = 178.0, 
92.0 X 10-^ 



4.7, 



(3.9) 



The (one-loop improved) mass of Higgs particles, charginos and neutralinos 
in the set (A) and (B) is given in Table 3 and 4, respectively. 

The set (A) is almost the same as the SPAla' parameter set, in which not 
only the heavier chargino xt but also the lighter chargino xt can decay into 
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tan/3 


M 


Ml 


Ma 


Ms 




10.00 


399.15 


100.13 


157.53 


610 


431 





771^2 




"ij. 








COS0U 


cos 64 


cos dg 


506.48 


524.14 


506.07 


530.14 


163.22 


187.37 


169.64 


9.4x10"^ 


8.5x10"* 


9.1x10"^ 


mei 






"IS2 








COS^c 


cos 0s 


cos 0f_i 


506.47 


524.16 


506.07 


530.14 


163.19 


187.38 


169.64 


0.033 


i.exiO"'' 


0.019 




"its 












cos^t 


COS Ob 


cos 0T 


345.37 


556.78 


469.43 


507.15 


150.07 


190.39 


170.02 


0.5567 


0.9266 


0.271 



Table 2: The value of the MSSM input parameters for set (B) 











107.12 


425.30 


424.90 


432.75 



xt 


xt 


X? 


xl 


X3 


xl 


184.2 


421.2 


97.75 


184.62 


398.30 


413.39 



Table 3: The pole mass of Higgs, charginos and neutralinos in set (A) 













122.50 


431.40 


431.0 


438.73 


xt 


xt 


X? 


X'2 


X^ 


xl 


147.08 


418.8 


97.61 


147.4 


404.0 


418.8 



Table 4: The pole mass of Higgs, charginos and neutralinos in set (B) 
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various two bodies because the hghter chargino is heavier than some sleptons. 
For the set (B), on the other hand, the mass of the hghter chargino m~+ is 
smaUer than the mass of all sfermions as well as the sum + m^o^. 
This means that the lighter chargino xt cannot decay into any two bodies 
and Xi has only three body decay modes, ffx^. 

4 Numerical results 

The one-loop electroweak corrections on various two-body decay widths of 
the lighter xf and heavier chargino xt for the parameter set (A) are shown 
in Table 5 and 6, respectively, where we do not display the decay modes 
with small branching fraction Br< 0.1%. While all decays in Table 5 are 
the electroweak processes, the process xt ~^ bti in Tabic 6 gets both the 
electroweak and QCD corrections through the loop contributions and the 
photon/gluon emissions. We define F = Fq -|- 5T, where Fq and 5T are the 
improved Born decay width and the one-loop correction, respectively. Note 
that the improved Born decay width is different from the Born width. We 
obtain Fq by replacing the tree-level masses by the one-loop renormalized 
pole-masses presented in Tables 3 and 4 in the tree amplitudes. 







To (GeV) 




F (GeV) 




5F/Fo 


Br 


xt 




3.91 X 10" 


-2 


3.78 X 10- 


-2 


-3.3% 


50.11% 


xt 


J^tifit 


1.33 X 10" 


-4 


1.19 X 10" 


-4 


-10.2% 


0.16% 


xt 


— > T'^i>r 


1.47 X 10" 


-2 


1.48 X 10" 


-2 


+0.1% 


19.58% 


xt 




1.06 X 10" 


-2 


1.07 X 10" 


-2 


+1.0% 


14.24% 


xt 




1.06 X 10" 


-2 


1.07 X 10" 


-2 


+ 1.0% 


14.22% 


xt 


^w+xti 


9.65 X 10" 


-4 


1.28 X 10" 


-3 


+32.3% 


1.69% 



Table 5: One-loop corrections on xt decay widths for set (A) 

In Table 7, we show the one-loop electroweak and QCD corrections to the 

various 3-body decay widths of xt ^hc parameter set (B), for which the 
two-body decays of xt ^-^^ kincmatically forbidden. The three-body decay 
widths of xt ^"^^ iiot shown, since, decaying dominantly into two bodies, it 
has extremely small three-body decay branching ratios. Note that for the 
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X2 



To (GeV) r (GeV) dV/V^ 



+ 1.8%(QCD) 



Br 



Xt Urf^ 1.54 X 10-1 1.48 X 10-1 -3.9% 4.20% 

xt ^f.fj't 1-36 X 10-1 1.46 X 10-1 +7.5% 4.13% 

Xt u^et 1.36 X 10-1 1.46 x lO-i +7.6% 4.14% 

xt T+i)^ 6.89 X 10-2 5.70 x 10-^ -17.3% 1.61% 

xt yu+z>^ 4.33 X 10-2 5.38 x 10-^ +24.2% 1.52% 

xt e+i>e 4.32 x 10-^ 5.37 x 10-^ +24.4% 1.52% 

xt W+x^ 1-93 X 10-1 2.07 x lO-i +7.0% 5.85% 

xt W+xt 8.66 X 10-1 9.93 x lO-i +14.6% 28.12% 

xt Zxi 7.53 X 10-1 8.56 x lO-i +13.7% 24.26% 

xt h^xt 5.97 X 10-1 6.07 x lO-i +1.75% 17.20% 

v~ 1 if -8.9%(ELWK) 

bti 2.82 x 10-1 2.57 x 10-1 <^ _ ' 7.43% 



Table 6: One-loop corrections on xt decay widths for set (A) 



decay modes involving quarks Xi ~^ QQXi^ the electroweak and the QCD 
corrections are separately given. 

We should note that the lighter chargino xt is the next hghtest SUSY- 
particle (NLSP) in the set (B). Since NLSP must be first produced by acceler- 
ators, it is important to study the production processes and the experimental 
signals of the chargino xf- calculate the full-one-loop electroweak cor- 
rection for the chargino pair production e~^e~ — > xtxi ^it the future linear 
colliders ( See also [261 [271 [28] ) . Fig-l the energy dependence of the total 
cross section is shown for the lighter chargino pair production at e^e" col- 
liders for the set (B) (and set (A)). The one- loop electroweak correction is 
of order of —10%. We note that the radiative correction of this order can be 
detectable at the proposed linear colhder. 

By combining the production cross section (Fig.l) and the decay branch- 
ing ratios (Table 7) for the set (B), we obtain the one- loop corrected cross 
sections for the direct experimental signals. Fig. 2 shows the energy depen- 
dence of the cross section for the two types of the chargino signals, e^e"+ 
missing energies and 4-jets + missing energies, at e^e" colliders for the set 
(B). 
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To (GeV) 


r (GeV) 




Br 


xt 




4.42 X 10-*^ 


4.84 X 10-6 


+9.4% 


20.18% 


xt 


f^'^'^t^Xi 


4.42 X 10-*^ 


4.84 X 10-6 


+9.4% 


20.18% 


xt 


r+z/^X? 


6.46 X lO"*^ 


7.22 X 10-6 


H 


hll.8% 


30.09% 


xt 


udxi 


3.35 X 10-*^ 


3.55 X 10-6 




\ -0.2%(ELWK) 
1 +6.3%(QCD) 


14.81% 


xt 


csxi 


3.33 X 10-6 


3.54 X 10-6 




1 -0.2%(ELWK) 
[ +6.3%(QCD) 


14.74% 




Table 7: One-loop corrections on xt 


decay widths for set (B) 



5 Comments 

The total cross section is a sum of the one-loop corrected cross section and 
the cross section of the hard photon radiation. The former is negative and 
large, while the latter is positive and large, originating mainly from the initial 
radiation. In our calculation for the hard photon radiation, we do not set an 
energy cut in the upper limit nor any angle cut. 

In the calculation of the one-loop correction of the decay widths the sit- 
uation is almost the same except that in the evaluation of the decay widths 
the contribution of the real photon (gluon) emission from both initial and 
final states is important. The correction proportional to the fermion mass 
m/ is expected to be large for the third generation r, t and b, and an addi- 
tional enhancement can emerge in the case of large tan /3 for r and b. Note 
that we set tan/5 = 10 for both sets (A) and (B). We see these effects in the 
numerical results presented in Tables 5, 6 and 7, where the branching ratios 
of the T modes are larger than e and /z modes. 

The radiative correction on the chargino pair production in the SPAla' 
scenario (set (A)) has been studied by the Wien group [26j. Unfortunately, as 
they adopt different treatment of the photon emission correction, the direct 
comparison of our result for the full electroweak correction with theirs is 
not possible. For the comparison we extract the "weak" correction Acr„eak 
defined by 

AcTweak = O'clwk " 0"BORN * SqeD — Charcf'- (5-1) 

We find that the result of GRACE/SUSY-loop is consistent with the previous 
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Figure 1: Total cross section for chargino pair production e^e^ — > XiXi fo'^ 
the set (A) and set (B). Solid line and dotted line denotes duoop and ctborn, 
respectively. 

result [26] . 

In the parameter setting (A), the produced chargino dominantly decays 
into two bodies, because the mass of the chargino m-+ (= 184.2GeV) is 
larger than the mass of the lighter charged sleptons, , and rric,^ as well as 
mw + (see Tables 1 and 3). Since the chargino cannot decay into any 
squark in two-body decay modes and BR(x]'' W^^Xi ~^ QQ'Xi) is negligibly 
small (see Table 5), we cannot expect the signals with the quark-jets from 
the Xi pair production. We find also that BR(£i ixf) = 1 and BR(z/£ —>■ 
I'iXi) = 1 in the parameter setting (A). This means that the most plausible 
experimental signal of the chargino pair production is the lepton pair plus 
the large missing energies. The precise measurement of the energy and the 
momentum of the r leptons is particularly important because r+r^ signal is 
the dominant mode in this case. 
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Figure 2: Total cross section for e'^e'^ XiXi {^^^eXTli^^^eXi) (a) and 
e^e~ — > XiXi ~^ {qq'Xi){'1Q'Xi) (b)- Solid line and dotted line denotes auoop 
and ctborn, respectively. 

In the parameter setting (B), the chargino xt with m-+ = 147.0GeV is 
lighter than in set (A) in which m-+ = 184.2GeV. We find the apparent 
shift of the production peak to the lower energy due to the different chargino 
mass value used in the set (A) and (B). As has been discussed in the previous 
section, two-body decays of xf are kinematically forbidden in the set (B). 
We can use both types of signals, the lepton pair plus large missing energies 
and the quark-jets plus large missing energies, for the chargino detection. 

We have developed a tool for the full automatic one-loop calculation of 
the MSSM processes, GRACE/SUSY-loop, which is characterized by the gauge 
symmetric and the on-shell renormalization scheme, and includes the various 
self-consistency check schemes. It is certainly a useful tool for the present 
and future precise analyses of the experimental data. 



Acknowledgement 

We would like to thank Y. Shimizu, T. Kaneko, G. Belanger and F. Boudjema 
for fruitful discussions. This work was partly supported by Japan Society 



18 



for the Promotion of Science under Grant-in- Aid for Scientific Research (B) 
(No.17340085). 



19 



Appendix A. Expressions of conterterms 

In this Appendix, we list the expression of the counterterms and the renor- 
mahzation constants in terms of two-point functions. 

The wavefunction renormahzation constants are expanded in the one-loop 
order as 

= l + l^Zx, (A.l) 

yi/2 _ {l + \5ZxY, X = Y , . 

- \'^6ZxY, X^Y ^^-^^ 

We use the following abbreviations; 

Cw = cos 9w, cl = COS0L, cr = COS0/J5 Ca = COS a, C/3 = cos/5, 
Sw = sm9w, SL = sm(f)L, s/j = sin0/j, SQ, = sina, = sin/3,(A.3) 

where 0l and (pn are the mixing angles which diagonalize the chargino mass 
matrix (see (2.10) of [TT]), and a is the mixing angle in the CP even Higgs 
sector, 
gauge sector 



6Zn= Re\JVU^)^r-^^-^^^^^^^^^^^l (A.5) 
6Z.g- Ke{yv^^[}}) + jp^i J. lA.bj 

Other counterterms appearing in the gauge sector are expressed in terms of 
bZw, SZb and 6Zg. For example, 

9 

SI 

9' 

Sw 

6cw 

Cw 

20 



3 

6Zg — -6Zw, 


(A.7) 




(A.8) 




(A.9) 


9 9 






(A.IO) 



fermion sector 



Sruf = -Re^imf) -mfReT/fimf), (A.ll) 
5Zf = ReJl^mj) + ReEf{mf) 

+ 2m f[ ReT.f{mj) + mfReEy{mf) ], (A. 12) 
SZf = ReJ:^{mf) - Re^firuf) 

+ 2mf[ ReEfinif) + mfReEy{mf) ]. (A.13) 

where the selfenergy function of Dirac fermion / is decomposed as 

E^(^) ^ E?(g2)l + (g2)75 + E^iq'U + E^(g^)^75. (A.14) 
tadpole terms 

5%^T';r. ^ = l,2, (A.15) 

where 5% is the tadpole counterterms. 
Higgs sector 



where 



5Ml 



= -i?eEAOAo(Ml) + MlRei:'j^oAo{M\), 



C(l,2) 
C(l,3) 
C(2,2) 
C(2,3) 

C(3,2) 



24 
-24 



c(^,J) +; 



.ReS^o^(M2) 



C ReE' 



Mz 



CSX 
Cc2p5Y 



sin 2/5(sin2/3 + sin2Q;)' 



Q 2 2 



2/2 2\ 



.2 „2 



3s^c^ 



44 + 2Sc«Cc,S^C/3, 



= -s, 



1) + 4 + 2SaCaSpCi3], 



(A.16) 



)1 



, (A.17) 
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C{3,3) 



-[4(2sl - 1) + 4 + 2SaCaSf3Cp], 



C(4,2) = + [sl(2c}-l) + sl + 2sacaspcp]. 



C(4,3) 



[sj(24 - 1) + 4 + 2SaCaS/3C(3]. 



(A.18) 



and 
5X 



i?e[n^^(0) + 2cwsw 



n^z(o) 



+ 



(4 - s2^)n^^(M2,) - ct^n^z(M|) 



(A.19) 



(5y = 



sin (a — /3)(5M4 + cos (a + l3)M2SZz 



+ 
+ 



9 



2Mw 
9 



cos{a + sin^ia - P)]STho 



2M^ 



w 



sm{a — (3) cos^(q; — (3)5T} 



(A.20) 



with 



5Zz = 2clr5Zg - Scl^SZw - sI^SZb. (A.21) 



The tadpole counterterms for the physical Higgs bosons are defined by 



5Tho 



cos a sm a 
— sin a cos a 



( 6% ] 



The counterterm of tan /? is then determined by 

S tan (3 = 

while the gauge-boson mass counterterms are given by 



1 tan P{6Zh, - SZh, - 2^ + 2^), 

2 Vi V2 



(A.22) 



(A.23) 



5Ml 



M^[2SZg - 36 Zw + cos^ (36 Zh, + sin^ P6Zh, 



-2cos^/3^-2sin^/3^] 

r2 ft- V I „„„2 njrry , „• 2 



M|[5Zz + cos' (36 Zh, + sin' 

6V2 
V2 



2cos^/3^-2sin^/?^], 



22 



(A.24) 
(A.25) 



Chargino sector 





f 








v 



5/1 



■-i2 ^2 I ^2 ^2 
,2 „2 „2 ^2 



.2 „2 
,2 ^2 



1 (244-ci-4)/X 2slcl/X 

1 /r)^2 „2 „2 „2 \ / V 0^2 „2 / 



SlSr{A22 - eLll22) - ClCr{Au - Uu) 
^2 _ „2 • 

SLSR{An - nil) - ClCr{A22 - eLll22) 



(A.26) 
(A.27) 

(A.28) 



where 



^11 



^22 



4 - 4, (A.29) 
2m^+[ReJ:f[{m^^) + m^+i?eEf/(m^+)] + i?eEfi(m^+), (A.30) 

, , ^ , . cos /3 , ^ sin /3 , ^ ^ 

+ (CLCi?M2 + ClSrMw-^ + SlCrMw-^)0Zxw 

+{clSr\/2Mw cos f3 + SLSRi^)^6Zfj^ 

+{slCrV2Mw sin P + slSrIj)^5Z^^ 

1 3 3 
+clSrV2Mw cos P{-6Zh, + — - — ) 

2 g vi 

+SLCRV2MwsmP{l5ZH, + ^ - (A.31) 
z g V2 

cos j3 sin j3 

+ {slSrM2 - SlCrMw—^ - ClSrMw-^)oZx^ 

+{-SLCR\f2Mw cos (5 + CLCRjJty^SZjj^ 

-SLCij"\/2MH^cos/9(^(5Zii-^ + ^ - — ) 
2 g vi 

-CLSRV2Mwsm(3{UzH, + ^ - (A.32) 

2 g V2 

2m\+[ReT.l[{m^+)+m^+ReT.Xl{m^+)] - i?eEfi(m^+),(A.33) 
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n 



m 



22 



X2 



[{sl + sl)5Z,. + cldZf,^ + ci^Z^J 



2 

—i?eS22 — m-+ReTior,(m 



V 

22V"'X^ 



(A.34) 



Neutralino sector 
1 



SZx 
6Mi 



[5Zn - 0^2^Zx^o - O^^SZg^ - O^^SZf^^] 



1 r6mii 



(A.35) 



Q2 



11 



'^l (P,9)7^(l,l) ip,q) 



(A.36) 



where (9 is the orthogonal matrix which diagonahzes the neutrahno mass 
matrix 

( nif. 



'"ni 



ruf. 



ruf 



(A.37) 



and 



6M, 



N 



( 5Mi -Mz svi/cos/5Ai3 Mz siysin/3Ai4 
* 5M2 Mz cw cos /?A23 —Mz cw sin /3A24 



* 



* 











(A.38) 



with 



Ai3 = 


5M| 

2M| ""^^ 


9 


S9\ 
9'' 


— tan [3 cos^ [35 tan /?, 


(A.39) 


Ai4 = 


5M| 

2M| ""^^ 




S9' 
9'' 


+ cot /3 cos^ [35 tan /?, 


(A.40) 


A23 = 


2M| + 




S9' 
9'' 


— tan [3 cos^ tan /3, 


(A.41) 


A24 = 


(5M| 2 
2M| + 


^9 
9 


S9' 
9'^ 


+ cot /3 cos^ [35 tan /3. 


(A.42) 



The phase factor rji in ( ]A.36ll is to convert the negative mass eigenvalue in 
( 1A.37I1 to positive, 

f 1 nin. > / . 
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sfermion sector 

For simplicity, we show only the expression for the first generation. 

6mj = -i?eSj^-(m^-), / = -Ui, ^2, c?2, ei, 62, z/g, (A.44) 

SZff = E'^jlmj), f = ui,U2,di,d2,ei,e2,iye, (A.45) 

o^^fifi ^ T' f = '^^d,e. (A.46) 

2 ^'^^ rrif — rrif 

Ji J 3 



SO, 



Smo^ - S{M^ cos 2/3 - m^) - cos^ OeSmj^ - sin^ 9e8m\^ 



sin 26*6 (ml^ - m|J 

(A.47) 

59^ = l ^"^"^(<) + ^';^^^«) (A.48) 
2 m4 - 



(5(cos^ + sin^ ^uW?^ - cos 2/3 - to^ + m^) - cos^ OdSm^ - sin^ Oddm^ 

~ '■ — ^^Tw 2 2~\ ' ~ 

sm20d{m~ — m- ) 

0,2 ai 



(A.49) 



QCD sector 



5Ms = -mg[ReEl{mg)+Re^^{m-g)], (A.50) 
= Rei:y{mg)+2ml[ReEl'{mg)+ReJ:y'{mg)], (A.51) 
<^^giuon = i?en;^(0), (A.52) 

^^^^ = -C[^-7E + M47r)]. (A.53) 
where C is the finite constant appearing at the one-loop vertex correction as 
i^V^ux)[C{j^ -1E + HiTT)) + •••]. (A.54) 
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Appendix B. External wavefunction renormal- 
ization constants 

We list the external wavefunction renormalization constant which ap- 

pears in the amplitude as 

M - -SZ'""' X (Born amplitude) (B.l) 

Gauge bosons 



Higgs sector 



Il'^{M^)=Il'^iM^)-6Zw, 
U'zziMl) = ll!zz{Ml) - bZzz. 



(B.2) 
(B.3) 



5Z^"o* = Sho//o(M|o). (B.4) 

Since the pole mass of /i" and in one-loop order does not agree with the 
tree mass, some complication appears. 



8Z^^ 



tth(Mlo,mlo,mlo) (B.5) 



where 



^hh 



_d_ 
9^ 



d 



dq^^q^-mjjo-^HH{q') 



(B.6) 



and M/jO is the one-loop improved pole mass of h^. The renormalized selfen- 
ergy functions are defined by 



Ihh 



^{q^)hh + 5Ml^ - q^{sin^ aReSZn, + cos^ aReSZn^), (B.7) 



E{q'^)HH = ^{q'^)HH + ^MIjj - q^{cos^ aRedZn, + sin^ aReSZH^), (B.8) 
t{q^)Hh = ^{q'')Hh + SMl^ - q"" cos a sin a{Re6ZH, - RedZn,), (B.9) 
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with 



cos2(a - 13)5 M\ 



(B.IO) 



g 



2Mw 
9 



cos(a — j3) sin^(a — I3)T^^°^ 



loop 



2Mw 

r2 \„-2i 



sm(a — + COS (a — l3)\T^o 

+M|[sm^(« + (3){5Zz + + 5Zh,) 

+ sin(Q; + (3) sin(a — — 5Zh2)\ 

+2 sin(a + /?) [sin (3 cos /? cos(q; + (3) — sin a cos /5]M| 



5vi 

Vl 



r2 



HH 
Hh 



—2 sin(a; + /5)[sin {3 cos /? cos(a + /?) + cos a sin /5]Mv , 

-ReT.HH{Ml) + M^o(cos2 a^Z^i + sin^ a^^/zj, 
— sin(/5 — a) cos(/5 — a)5M\ 



(B.ll) 
(B.12) 



- sin3(/3 - «)tS°p + ^ cos3(a - P)T^,^^ 



2M, 



-M| [sin (a + (3) cos(a + + 5Zh, + ^Z^J 

+ sin a cos a{5Z Hi — SZh.^)] 

A . sin 2a cos 2q; sin 2/3, , ^9 
+|_„„(2„ + 2,5)(1 + 00= 2,3) _p_^]M|- 

.1 . sin 2a COS 2a sin 2/9 2 "^^2 



+|-=i„(2« + 2/3)(l-coB2^)+ ^j„p„-2^) 
The expressions (IB.SP and (]B.6P agree with those given in [25] . 



^"2 



where 



5 ' , 2 

— 2-.H±H±[(l 



dq 



Af2 , 



[B.13) 



[B.14) 



d 



The renormahzed selfenergy functions appearing in flB.14p are given by 
Hq^)h±h± = W)h±h± + 5Mi±^± - q^slReSZn, + clRedZn,), (B.15) 
S(g')//±G± = S(g2)^,±G± + 5M^±G± - q^Cf3Sp{Re5ZH, - ReSZnJ, (B.16) 
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with 



SMjj±H± = SMloAo+M^{6Z,-2cf-^-2sf-^), (B.17) 
ml^a^ = 5M2o^o-c^s^M^(^-^). (B.18) 



In the one-loop order, (IB.SP and (]B.14p become 



6Zlt = KoAM^o), (B.19) 
= S'^i^4Mi±). (B.20) 



Chargino 



= 2M,+ (Sf (M|,) + M,+ Sr(M|,)) + Er(M|,) - -{6Zg + SZg), 

(B.21) 

where the chargino selfenergy functions are decomposed as (]A.14I) . In terms 
of the renormalization constants introduced in section 2, the chargino wave- 
function renormalization constants are given by 

5Z22 = smcp'^dZx^ + cos(l)\5Zjj^, 

5Zg = sin (pl5 Z + cos (plSZji^. (B.22) 

Neutralino 

5Z5f = 2M^o[Ef/(M|o) + M^oE]^'(M|o)] + S,^(M|o) - SZ^. (B.23) 

where i = 2,3,4 and the selfenergy function of Majorana particles is decom- 
posed as 

j:M^^fifn + ^^i9')l (B.24) 

The neutralino wavefunction renormalization constant appearing in (]B.23P 
is expressed in terms of the basic renormalization constants introduced in 
section 2 as 

SZii = ^{ReSZ)k{ON)ik{ON)ik, no sum over i, (B.25) 

k 

where 

SZ, = iSZx,5Zx^,SZfj^,SZfj^). (B.26) 
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and O is the orthogonal matrix which diagonahzes the neutrahno mass ma- 
trix. See (IX37D . 

Note that for unstable particles, even if the residue condition is imposed 
on the propagator at the pole position, there is a non-vanishing SZ^^^ which 
is ultraviolet-finite and purely imaginary. For example, we can easily check 
that which is obtained from (lB.2ip by changing the index 2 to 1, is 

purely imaginary. We can neglect such contributions if perturbation works. 

References 

[1] For review see, for example, P. Fayet and S. Ferrara, Phys. Rep. 32 
(1977) 249; H. Haber and G. Kane, Phys. Rep. 117 (1985) 75. 

[2] H.P. Nilles, Phys. Rep. 110 (1984) 1. 

[3] S.L. Glashow, Nucl. Phys. 22 (1961) 579; S. Weinberg, Phys. Rev. Lett. 
19 (1967) 1264; A. Salam, Proceedings of the Nobel symposium, 1968, 
Lerum, Sweden. 

[4] J. Fujimoto et al, Comp. Phys. Comm. Ill (1998) 185. 

[5] J. Fujimoto et al, Comp. Phys. Comm. 153 (2003) 106. 

[6] J. A. Aguilar-Saavedra et al, Eur. Phys. J. C46 (2006) 43. 

[7] T. Ishikawa et al, KEK Report 92-19, 1993, The GRACE manual 
vl.O; F. Yuasa al, Prog. Theor. Phys. Suppl. 138 (2000) 18, 
hep-ph/ 0007053|. 

[8] J. Fujimoto, T. Ishikawa, M. Jimbo, T. Kaneko, T. Kon and M. Kuroda, 
Proceedings of the XVIIth International Workshop on High Energy 
Physics and Quantum Field Theory, (2004) 26, | hep-ph/0402l44| . 

[9] J. Fujimoto, T. Ishikawa, M. Jimbo, T. Kon and M. Kuroda, Nucl. 
lustrum. Methods Phys. Res. A 534 (2004) 246. 

[10] J. Fujimoto et al, Nucl. Phys. Proc. Suppl. 157 (2006) 157. 

[11] M. Kuroda, KEK-CP 080 (1999), |hep-ph /9902340|. 



29 



[12] K. Hikasa, SUSY manuscript version July 5, 1995, (1995), [unpublished]. 

[13] J. Rosiek, Phys. Rev. D41 (1990) 3464; erratum KA-TA-8-1995. 

[14] P.H. Chankowski, S. Pokorski and J. Rosiek, Nucl. Phys. B423 (1994) 
437. 

[15] Y. Yamada, Phys. Lett B530 (2002) 174. 

[16] J. Guasch, J. Sola and W. Hollik, Phys. Lett. B437 (1998) 88. 

[17] G. Belanger et ai, Phys. Rep. 430 (2006) 117. 

[18] A. Yamada, Phys. Lett. B263 (1991) 233; Z. Phys. C61 (1994) 247. 

[19] A. Dabelstein, Z. Phys. C67 (1995) 495. 

[20] D. Pierce and A. Papadopoulos, Phys. Rev. D50 (1994) 565; Nucl. Phys. 
B430 (1994) 278. 

[21] S. Alam, K. Hagiwara, S. Kanemura, R. Szalapski and Y. Umeda, Phys. 
Rev. D62 (2000) 095011. 

[22] W. Hollik, E. Kraus and M. Roth, C. Rupp, K. Sibold, D. Stockinger, 
Nucl. Phys. B639 (2002) 3. 

[23] J. Guasch, W. Hollik, J. Sola, JHEP 0210 (2002) 040, [hep-ph/0207364] . 

[24] T. Fritzsche and W. Hollik, Eur. Phys. J. C24 (2002) 619. 

[25] S. Heinemeyer, W. Hollik, J. Rosiek and G. Weiglein, Eur. Phys. J. C19 
(2001) 535. 

[26] W. Oiler, H. Eberl and W. Majerotto, Phys. Rev. D71 (2005) 115002. 

[27] W. Kilhan, J. Renter and T. Robens, Eur. Phys. J. C48 (2006) 389. 

[28] W. Kihan, J. Renter and T. Robens, DESY 06-189 (2006), 
|hep-ph/0610425| . 



30 



